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Introduction 


“Oil sands” is a type of petroleum deposit in which the oil is very 
thick and sticky (called “bitumen”) and mixed with sand/water/clay 


Alberta’s oil sands areas 


Province 
of Alberta 





- Bitumen found close to the surface 
may be mined; deeper deposits need 
to be heated and then pumped to 
surface 


i surface mining region 


from Energy Resources Conservation Board, 2011 



Proven reserve of ~170 billion barrels 

Production expected to double by 2020 

Additional monitoring needed to better understand the emissions of 
the oil sands region and its impacts 

- Joint Canada and Alberta plan for monitoring of the air, water, and wildlife in 
and around the oil sands 

- Satellites provide large scale spatial and temporal coverage and extent 
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What insights can Aura provide on Air Quality 

in the Oil Sands? 


Aura-OMI 

(Ozone Monitoring Instrument) 

OMI provides tropospheric vertical 

column densities (VCDs) 

• N0 2 : NASA SP v2.1 & KNMI 
DOMINO v2.0 

• S0 2 : NASA PCA*-beta release 

Largely eliminates artifacts and 2x reduction 
in noise compared to operational product 

• For all products original AMFs 
replaced by new Environment Canada 
AMFs that are based on higher 
resolution input data [McLinden et al., 
ACP, 2014] 

• SP and DOMINO datasets combined 
since remaining difference is primarily 
stratospheric N0 2 removal, and 
combined data appears to work best 
over Canada 

• Principle Component Analysis method, Li et al., GRL, 

2013, product being evaluated; See Nick Krotkov talk 


Aura-TES 

(Tropospheric Emissions Spectrometer) 


TES provides a volume mixing ratio 
(VMR) profile 

• Used recent Version 6 Lite products* 

- New CH 3 OH and HCOOH products 

- significant amount has been reprocessed 

• NH 3 , CH 3 OH, HCOOH 

- Peak sensitivity varies between 1 -2 km 

- Typically 1 DOFS or less 

■ Not much vertical “profiling” 

■ Reported as a RVMR 

- Boundary layer weighted averaged VMR 
value where TES is most sensitive 

• CO 

- Peak sensitivity typically ~3-km 

- Typically 1-2 DOFS 

- For comparison purposes we report the 
VMR at the peak vertical sensitivity in the 
troposphere defined by the averaging kernel 
(AK Peak) 
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OMI view of Central Canada 
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Evolution 
of OMI NO, 



Movie goes here 
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Evolution 
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2005-2007 


2011-2013 



2005 - 2013 E(NPRI*) = 100 kt[S0 2 ]/yr 

average *NPRI = National Pollutant Release Inventory 
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OMI S0 2 over the oil sands 


75 kt 
— 22 kt 


VCD [DU] 


S0 2 emissions due to 
upgrading 

• converting bitumen to 
synthetic crude 

Only two significant S0 2 
point sources, both in 
southern [S] mining region 
Northern [N] mines pipe 
bitumen off-site for upgrading 


E(OMI) = 89 kt[S0 2 ]/yr 


3-year averages 
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Emissions and lifetime 
determined by fitting 
the downwind decay of 
S0 2 

[similar to Beirle et al.. Science, 
2011] 


2011 2012 2013 

NPRI = National Pollutant Release Inventory 









2005-2007 


2011-2013 


Mining 
Regions 



OMI N0 2 over the oil sands 


N0 2 emissions sources : 

• upgrading (50%) 

• large vehicles (50%, more uncertain) 

Two significant point sources (upgraders), 
both in southern [S] mining region. 
Significant area sources in [S] and [N] 
Change in distribution consistent with 
expansion into the [N] 


VCD [10 15 cm 2 ] 


2005 - 2013 

E(NPRI*) = 53 kt[NO z ]/yr 

E(OMI) = 55 kt[NO z ]/yr 

average 

*NPRI = National Pollutant Release Inventory 



Emissions and lifetime determined by 
fitting the downwind decay of NO z 

[similar to Beirle et al., Science, 2011] 

NOx / N0 2 = 1.35 from AQ model 


Is this difference real, from sources not 
reported to the NPRI (e.g., construction)? 


2006 2007 2008 2009 2010 2011 2012 2013 NPRI = National Pollutant Release Inventory 

Year 









TES Infrared Satellite Observations 


TES Global Survey (GS) Mode 

- Nadir pointing (16-day repeat cycle) 

- Spacing of ~180 km along track 

- No Global surveys taken after 2011 



TES Special Observation (SO) Mode 
* Higher sampling density over shorter tracks 
- Transect: regional pollution studies 
- Over the Oil Sands: 

- Begun July 14, 2012 (over 2-years) 

- Observations every 2-7 days 

» Over 125 SO to date 

- 20 contiguous targets 

- 12-km sampling along track 

- Between 56-58°N covering 240 km 

» Centred on the oil sands 

- Each target is 5x8 km 



Used recently produced TES Version 6 Lite Products (Susan Kulawik). 
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Ammonia : Central Canada 


TES NH3 SeasomAII (ppbv) 



- Period from 2004-2014 


■ mostly GS before 201 1 and SO over the OS after 

- 2x2° grid averages 

- Overall spatial gradient consistent with NPRI emissions database 

- Potential decrease in NH 3 over the oil sands region? 
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CO HCOOH CH30H NH3 

VMR [AK Peak] (ppmv) RTVMR (ppbv) RTVMR (ppbv) RTVMR (ppbv) 


Trends in TES Over Oil Sands Region? 


TES NH3 Season:AII Lat= 56:58 Lon=-110:-112 



TES CH30H Season:AII Lat=56:58 Lon=-110:-112 
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Special obs (SO), 2012-2014 

# Surface mining 


Approach: compare global survey 
(2004-2011) with special observations 
(2012-2014) over oil sands region 
Challenging: Not many GS values over 
oil sands region (2004-2011) 

Any trends? 

• Large values in summer 2012 
potentially due to biomass burning? 

• Need to identify (and filter) 
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CO HCOOH CH30H NH3 

VMR [AK Peak] (ppmv) RTVMR (ppbv) RTVMR (ppbv) RTVMR (ppbv) 


Are OS Values Different than Across Central Canada? 


TES NH3 Season:AII Lat=54:60 Lon=-118:-98 
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Approach: compare global survey oil 

sands with larger central Canada region 

• Seasonal cycle in GS and SO 

• Global survey (2004-2011 ) values 
used to define “typical” values from a 
large region in central Canada. 

• GS (2004-2011 ) values over central 
Canada are similar to the (relatively 
few) GS values over OS 
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CO HCOOH CH30H NH3 

VMR [AK Peak] (ppmv) RTVMR (ppbv) RTVMR (ppbv) RTVMR (ppbv) 


Are OS Values Different than Across Central Canada? 


TES NH3 Season:AII Lat=54:60 Lon=-118:-98 
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Approach: compare global survey oil 

sands with larger central Canada region 

• Seasonal cycle in GS and SO 

• Global survey (2004-2011 ) values 
used to define “typical” values from a 
large region in central Canada. 

• GS (2004-2011 ) values over central 
Canada are similar to the (relatively 
few) GS values over OS 
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“Climatological” Oil Sand Region Satellite Observations: 

Period 2012/07/14 - 2014/04/03 (Warm Season: M-J-J-A-S) 

Ammonia (NH 3 ) Methanol (CH 3 OH) 
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Satellite Validation: 2013 Intensive Oil Sands Field Campaign 


Sept. 3 rd , 2013 : Flight 18 


Sept. 5 th , 201 3: Flight 20 
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* Dedicated TES overpass spirals 

* Clear conditions 


Canada 






Satellite Validation: 

2013 Intensive Oil Sands Field Campaign 


* Comparisons of TES and OIV 

- Period from Aug. to Sept. 2012 

- Aircraft 

- Surface 

■ In-situ / Remote (Pandora) 

- AQ model 

■ GEM-MACH: 2.5 x 2.5 km 

• Validation of new TES CH 3 OI 
and HCOOH products 


• Very preliminary TES/aircrafl 
comparison results show: 

- -20% NH 3 

- -30% CH 3 OH 

- -40% HCOOH 

- -10% CO 

* Waiting on QC for aircraft 
N0 2 and S0 2 observations 
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Profiles Differences 
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Final Remarks 


* The 10-years of the Aura OMI and TES satellite observations are 
providing valuable insight on the air quality in and around the 
Canadian oil sands region. 

* Some highlights presented include: 

- OMI sees clear enhancements in N0 2 and S0 2 over the oil sands 

■ comparable with medium-sized city (~1 M) or large power plant 

■ distributions are consistent with location of sources 

- N0 2 increasing 

■ possibly at a rate faster than NPRI emissions would suggest 

- S0 2 showing slight decline and consistent with NPRI 

- S0 2 analysis : possible due to improvements due to new PCA algorithm 

- Initial analysis of NH 3 , CO, and VOCs (CH 3 OH, and HCOOH) indicates: 

■ TES does not detect large elevated concentrations directly over the oil sands 
mining regions 

■ Potentially a decrease of NH 3 over oil sands region 

- In the presence of sulphur and NOx : NH 3 (gas) — ► NH 4 + (aerosol) 

- Initial TES/aircraft validations show general agreement 
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Status and Availability of Infrared Satellite Obs 

TES Special Observations over the OS for the past ~2 years 


Over -130 special 
observations over 
the oil sands to 
date 

- Measurement 
every 2- 7 days 


OS Field Study 
• 9 overpasses 

TES transects 
of oil sands 
begun July 14, 
2012 
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NPRI ^-National Pollutant Release Inventory 

OMI 


3-year averages 



Trend = 3.3+ 0.4%/yr 
Trend = 5.7 ± 0.8%/yr 


E and x determined by fitting the 
downwind decay of N0 2 [similar to 
Beirle et al.. Science. 20111 

NOx / N0 2 = 1.35 assumed 


Is this difference real, from sources not 
reported to the NPRI (e.g., construction)? 
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• Effective lifetime short, reflecting the 
very rapid drop-off in N02 from its 
source (near background ~40 km 
away) - suggests higher OH levels 

• Some evidence for an increase in 
lifetime as N0 2 increases? This would 
suggest a transition towards a VOC- 
limited regime 
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Nice Features of the OE Approach 


Retrieval Errors 

- Straight-forward method of estimating 
retrievals errors 

- E={K T S m l K+S a - l }~ 1 ^ total error 

Averaging Kernels (A) 

- A^KTSJK + Sa 1 }- 1 KJSJK 

- Describes the relative weighting of the 
retrieved product, x r to the “true” atmosphere, 

x, and a priori x a 

x r = Ax + (I - A) x a 



Ideal Retrieval 
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DOFS (degrees-of-freedom for signal): 

Number of independent pieces of information in 

the measurement. 

DOFS=trace(A) 

Estimate of the vertical resolution: FWHM 

AK varies from profile-to-profile depending: 

Instrument (i.e. noise, nadir/limb viewing) 

Atmospheric state 

- i.e. temperature, trace gases, clouds 
Constraints 
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Comparison Methods : 

Apply the Observational Operator 

Provides the best “apples-to-apples comparisons 

- Requires a comparison profile 

Comparisons the satellite measurement information only 

- Essentially how TES would “see” the atmosphere measured by the aircraft 

- Put the high resolution data (aircraft) on the low resolution (satellite) 

- Removes the influence of the retrieval a priori when subtracted 


Satellite . x r A. sa ^x ( I A. sa ^)x a 


Air (TES obs) . X a i r ^sat^air "J" ^sat)^a 



Profiles 


Differences 
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Comparison Methods: 

Representative Volume Mixing Ratio (RVMR) 


* Collapse all information to a subset of 
level(s) where the retrieval is most 
sensitive 


• Reduces the influence of the a priori 

* Useful for retrievals with limited 
information 






• ~1 DOFS or less 

Useful generating maps, or comparing with 
non-profile single values (i.e. surface) 

Can be thought of as a “boundary layer” 
weighted average VMR where the satellite 
is most sensitive. 

As an example : 

• TES is most sensitive to NH 3 ~800 mb 

• DOFS = 0.83 

• RVMR = ~5.0 ppbv 

• ~2-3 km vertical resolution 

• Note: little sensitivity at surface 
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Challenges for Minor Trace Species: 

Example Ammonia (NH 3 ) from TES 



Relatively Weak Atmospheric IR 
Signal 

~ tenth’s to a couple degrees BT signal 
compared to a background of -300 K 

Detectability 

• - 0.5 - 1 ppbv under ideal conditions 

• thermal contrast plays a role 


NH Averaging Kernel NH (ppb) 

• TES is most sensitive to NH 3 between 900 
and 700 mb 

• ~2 km vertical resolution 

• 1 piece of information or less: DOFS<1.0 

• Representative Volume Mixing Ratio (RVMR) 

• Collapse all information to a single point that 
represents the information content: 

• Easier to compare with in situ measurements 
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